A Generalization of the Preceding Paper U A Rabi 
Oscillation in Four and Five Level Systems" 

Kazuyuki FUJII * 
Department of Mathematical Sciences 
Yokohama City University 
Yokohama, 236-0027 
Japan 

Abstract 

In the preceding paper |quant-ph/ 03 12060 we considered a general model of an 
atom with n energy levels interacting with n-1 external laser fields and constructed 
a Rabi oscillation in the case of n =3, 4 and 5. 

In the paper we present a systematic method getting along with computer to 
construct a Rabi oscillation in the general case. 
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In the preceding paper jT] (and [2J) we considered a general model of an atom with n 
energy levels interacting with n-1 external laser fields and constructed a Rabi oscillation 
in the case of n = 3, 4 and 5. Concerning more realistic model on two level system see |3j 
and 

The purpose of this paper is to present a systematic method getting along with com- 
puter to construct a Rabi oscillation in the general case. 

To begin with, let us make a brief review of the model. We consider an atom with n 
energy levels {{\k), E^) | < k < n — 1} which interacts with n — 1 external fields. We 
set Afc = Ek — E for 1 < k < n — 1 and assume the condition 

Ei — E > E 2 — Ei > ■ ■ ■ > E n -i — -E„_2 

for simplicity. 

We subject the atom to n — 1 laser fields having the frequencies u>k equal to the energy 
differences — A^._x = Ek — E^-i- As an image see the following picture : 
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Then the evolution operator U (t) defined by the Schrodinger equation 

i^-U(t) = HU(t) (h = l), 
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where H is the Hamiltonian given in [T| (we don't repeat it here), is given by 



U(t) = e - ltE °V^e- ltc , 



[I) 



where V = V(t) is 
/ 
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and the constant matrix C consisting of all coupling constants 



/ 



C= C(gx,g 2 ,---,9' 
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Therefore the remaining problem is to calculate the exponential e~ ttc , which is however 
very hard. 

In pQ we determined it by use of diagonalization method of matrices in the case of n 
= 3, 4 and 5. If all coupling constants are equal {g\ = g 2 = • • • = 9W-i) then the situation 
becomes very easy, see for example 0. In the following we present a method to calculate 
e -itc j n g enera i case, which is universal in a sense. 



Now let us review the formula in [H] within our necessity. Let A be a matrix in 
M(n, C) and {Ai, A2, • ■ ■ , A n } the set of eigenvalues (containing multiplicities) of A. Then 
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the exponential e ltA is given by the clear formula 

e~ itA = f (t)l n + f\(t)A + f 2 (t)A 2 + ■■■ + fn-^A"' 1 (3) 

with 

fi(t) = (-l) n+1 £ ^ n - 1 - 1 }* for 0</<n-l. (4) 

Here {pi,p 2 , • • • , p n -i,Pn} are a kind of fundamental symmetric polynomials of {Ai, A 2 , • • • , A„} 
and {(pi)fc, (p 2 ) k , • • • , (Pn-i)*} are them consisting of {Ai, • • • , A fc _i, A fc+ i, • • • , A n }. For ex- 
ample, 
n=3 : 

f,{t)l, + h{t)A + f 2 {t)A 2 (5) 



with 



n=4 : 



with 



e -itA = 



f ( t ) = A 2 A 3 e- itAl AiA 3 e- itA2 AiA 2 e~ ifA3 

/oU (A 2 -A 1 )(A 3 -A 1 ) (A 1 -A 2 )(A 3 -A 2 ) (Ai - A 3 )(A 2 - A 3 ) ' 

f() (A 2 + A 3 )e- aAl (Ax + A 3 )e- ftAa (Ai + A 2 )e~^ 

M J " (A 2 - Ai)(A 3 - Ax) (Ax - A 2 )(A 3 - A 2 ) (Ax - A 3 )(A 2 - A 3 ) 

hit) = (A 2 -Ax)(A 3 -Ax) + (Ai-A 2 )(A 3 -A 2 ) + (A x - A 3 )(A 2 - A 3 ) ' 



e~ UA = f (t)U + hit) A + f 2 (t)A 2 + / 3 (t)A 3 (6) 



A 2 A 3 A 4 e- iai A 1 A 3 A 4 e- i ' A2 
/o(t) ~ (A 2 - Ai)(A 3 - Ai)(A 4 - Ai) + (Ai- A 2 )(A 3 - A 2 )(A 4 - A 2 ) 
A 1 A 2 A 4 e- itA3 AAaAge"^ 4 
+ (Ai - A 3 )(A 2 - A 3 )(A 4 - A 3 ) + (Ax - A 4 )(A 2 - A 4 )(A 3 - A 4 )' 
() (A 2 A 3 + A 2 A 4 + A 3 A 4 )e-^ (AiA 3 + AjA 4 + A 3 A 4 )e~ aA2 
M j " (A 2 - A0(A 3 - A0(A 4 - AO (Ax - A 2 )(A 3 - A 2 )(A 4 - A 2 ) 
(AjA 2 + AxA 4 + A 2 A 4 )e~^ 3 (A X A 2 + AxA 3 + A 2 A 3 ) e -^ 
( Ax - A 3 ) ( A 2 - A 3 ) (A 4 - A 3 ) (Ax - A 4 ) ( A 2 - A 4 ) ( A 3 - A 4 ) ' 
, ( ) = (A 2 + A 3 + A 4 )e~^ (Ax + A 3 + A 4 )e~^ 

M } (A 2 - Ax)(A 3 - Ax)(A 4 - Ax) + (A x - A 2 )(A 3 - A 2 )(A 4 - A 2 ) 
(Ax + A 2 + A 4 )e-^ (Ax + A 2 + X 3 )e' u ^ 

(Ax - A 3 )(A 2 - A 3 )(A 4 - A 3 ) + (Ax - A 4 )(A 2 - A 4 )(A 3 - A 4 )' 
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J o\"J 



(A 2 - A X )(A 3 - Ai)(A4 - A x ) (Ai - A 2 )(A 3 - A 2 )(A 4 - A 2 ) 



-it\ 3 



-UX4 



( Ai - A 3 ) ( A 2 - A 3 ) ( A 4 - A 3 ) (Ai - A 4 )(A 2 — A 4 )(A 3 - A 4 ) 

It is notable that our formula (J3J) with (jlj) is convenient to use computer if the eigen- 
values are known. 



Next, what we do is to look for the eigenvalues of C in (j2j) in order to use the formula 
above. The characteristic polynomial of C is 

A -g\ 

~9i A -g 2 

-92 A -g 3 

f n (\) = \\l n -C\= ... (7) 

—gn-2 A —g n -i 
-g n -\ A 

By using the Laplace expansion of determinant it is easy to see 

/n(A) = A/ n _ x (A) - gUfn-2(X); / (A) = 1, /i(A) = A. (8 

For example, 
f 2 (\) = \ 2 - g 2 , 

/ 3 (a) = a(a 2 -< ?1 2 -< 72 2 ), 

/ 4 (A) = A 4 -^ 2 + ^ 2 + ^ 2 )A 2 + ^ 2 , 

/ 5 (A) = A (A 4 - {g 2 + g\ + g 2 + £ 2 )A 2 + {g 2 g 2 3 + g\g\ + g 2 2 gl)) , 

/ 6 (A) = A 6 - (g\ + g\ + g\ + g\ + g 2 5 )X 4 + (g 2 g 2 3 + g\g\ + g\gl + g\g\ + g\g\ + g 2 3 g!)X 2 

222 
-919395 

A (A) = A (A 6 - (gf + <? 2 + gl + <? 4 2 + # 5 2 + <? 6 2 )A 4 + + g\g\ + £ 2 <7 5 2 + <? 2 # 6 2 + g 2 2 gj+ 

22, 26, 22, 22, 2 2\\2 / 2 2 2 , 222, 222, 22 2\\ 

#2#5 + 929 6 + 9 3 9 5 + 9s96 + 949 6 )^ ~ (#i£ 3 #5 + 9i9 3 9& + 9i9±9% + 929i9&) ■ 



(9) 



Now, let us look for the general form of the characteristic polynomial of /„(A). 

Main Result 

(i) n = 2m : 

/ 2m (A) = A 2m - 2 A 2m " 2 + • • • + (-1) k fo k \ 2m - 2k + ■■■ + (-l) m ~ V 2m - 2 A 2 + (-l)"> 2m (10) 
with 

2m- 1 
02 = Yl 9i, 



i=l 

2m-l 



, v— v 2 2 2 2 

k* = E 9i 1 9i 2 ---9i k _ 1 9 lk 

l<ii<i2<h<-- ■ <ik-i<ik\ |*l — »2|>2,|i2-«3|>2,---,|ifc_i-ifc|>2 

for 2<fc<m-l. (11) 



fii) n = 2m + 1 : 



iWi(A) = A {A 2m - cp 2 X 2m ~ 2 + ■■■ + (-l) k cp 2k X 2m - 2k + ■■■ + (-1)— V 2 m- 2 A 2 + (-l)> 2m } 

(12) 



with 



2m 



V?2 = Y,9 2 i, 



i=l 

2m 



^2k= E ahl-'-at-A 

l<ii<i2<«3<— <ifc-i<»fc; |»i-i2|>2,|i2— »3|>2,— ,|ifc_i— ifc|>2 

for 2 < k < m. (13) 



The key point of the proof is the following equations. 

2m- 1 

l<ii<i'2<ki<---<ik-i<ik; »2|>2,|i2-i3|>2,— ,|ifc_i-ifc|>2 
2m-3 

E 2 2 2 2 | 

9i 1 9i 2 ' ' ' 9i k _ 1 92m-l ' 

l<ii<j 2 <i3<---<ifc-i; I»i-i2|>2,|i2— »3|>2,— ,|ifc_2— »fc-i|>2 
2m-2 

E 2 2 2 2 

9i x 9i 2 9i k _ 1 9i k 

l<ii<i2<h<---<ik-i<ik\ »2|>2,|z2-i3|>2,— ,|ifc_i-ifc|>2 
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and 



2m 

V2k= E g\g\- ■ ■ g\^g\ 

l<ji<i 2 <«3<---<*fc-i<«fc; \h— i2]>2,|i2— *3|>2,— ifc[>2 
2m-2 

E2 2 2 2 I 

fi^Sta ' ' ' 9i k _ 1 92m i 

l<ii<i2<*3<— <ifc-i; i2|>2,|ia— ts|>2,— ,|ffc_2— i*_i|>2 
2m- 1 

E 2 2 2 2 

fin 5^2 ' ' ' 9i k _ 1 9i k ■ 

l<ii<i2<«3< •••<ifc_i <ifc; *2|>2,|i2— i3|>2,—,|ijt_i— ifc|>2 

For the proof we have only to use the mathematical induction. We leave the remaining 
part to readers. 

A comment is in order. Since the matrix C in (J2J) is real symmetric (of course hermi- 
tian) its eigenvalues are all real. 

Next, let us solve the characteristic polynomial of C. From 
n = 2 : 

Ai = gi, A 2 = -g%. 

n = 3 : 



Ai = \/fii+fi 2 , A 2 = 0, \ 3 = -\/9i+ gi 

n = 4 : 

v^+v^B x/A-y/B VA-VB VJ+y/B 

Ai- - , A 2 - - , A 3 - - , A 4 - - 

where 

A = g\ + (gi + g 3 ) 2 , B = g 2 2 + ( 9l - g 3 ) 2 . 

n = 5 : 



+ y/A-y/B . n . yfA-y/B . v 7 ^ + 
Ai - , A 2 - , A 3 - U, A 4 , A 5 

where 



A = gl+gl+gl+gl+2\Jglgl + g\g\ + ghl B = gl+gl+gj+gj-2\J 'g\gl + g\g\ + g\g\. 
n = 6 : 
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By setting A 2 = x we can use the Cardano formula (see for example j7j). The solutions 
of the equation 

x 3 — ax 2 + bx — c = 

with 

2 i 2 i 2 i 2 i 2 i. 2 2 i 22, 22, 22, 22, 22 222 

a = g x + g 2 + 9 3 + 9 4 + g 5 , o = g x g 3 + g x g 4 + g x g 5 + g 2 g 4 + g 2 g 5 + g 3 g 5 , c = g±g 3 g 5 
are given by 

Xl = u + v + ^, x 2 = au + cr 2 v + |, x 3 = a 2 u + av + |, 
where a = e 2m ^ 3 and Uq, Vq are each solution of the binomial equations 



, — q + y/q 2 + 4p 3 , — g — V? 2 + 4p 3 
2 ' 2 

with 

b a 2 ab 2a 3 

P ~3~~9' q ~~ C+ Y~l>7' 

It is not difficult to show that we can choose x\ > x 2 > x 3 > 0. Therefore the solutions 
are 

Al = y/Xl, A2 = \fx 2 ^ A3 = \fX~3~, A4 = — y/xi, A5 = — y/x~2, Xq = —y/Xi. 

n = 7 : 

In the case of n = 6 we have only to change a, b and c to 

« = # 2 + g\ + g\ + g\ + g\ + g%, 

1 2 2 1 2 2 1 2 2 1 2 2 1 2 2 1 2 2 1 2 2 1 2 2 1 2 2 1 22 
= g x g 3 + + g x g 5 + g\g& + g 2 g^ + g 2 g^ + 9 2 g& + 9 3 9 5 + g 3 9 e + g^g 6 , 

222, 222, 222, 222 
C = #l#3#5 + £l#3#6 + 9ig^6 + 929496 

and obtain 

Al = \fxi, A2 = \f%2-) A3 = WXs, A4 = 0, A5 = —y/Xs, Xq = —\fx~2-, Xj = —\JX\. 



For n = 8 and 9 we can in principle solve the characteristic polynomial ((7J) by use of 
the Ferrari or Euler formula (see [Zj), which is of course very complicated. However, for 
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n > 10 it is impossible to obtain algebraic solutions of the characteristic polynomial by 
the famous Galois theory. Therefore we must appeal to some approximation method like 
the Newton's one (which is well-known), see the following picture. 

y = f( x ) 




Here for a given, the points a±, 02, etc are given by 

f(a) /(01) JK-i) 

di — a — , , <22 — at — 777 — -, • • • , ai- — ak-i — 777 r, • • • . 

/ (a) / (ai) / (a fc _i) 

For an appropriate number n we have only to set A = a n . By changing a in the general 
case we obtain a set of approximate solutions {Ai, A 2 , • • • , A n }. 



In this paper we generalized the result in pQ by making use of . If we can find exact 
or approximate solutions of the characteristic polynomial then we have the (exact or 
approximate) evolution operator by Q and (Jlj). That is, this means that a complicated 
unitary matrix in qudit theory was obtained. It may be possible to replace long quantum 
logic gates in qudit theory with few unitary matrices constructed in the paper, which will 
be discussed in another paper. 



References 

[1] K. Fujii, K. Higashida, R. Kato and Y. Wada : A Rabi Oscillation in Four and Five 



Level Systems, to appear in Yokohama Mathematical Journal, quant-ph/0312060 



9 



[2] K. Fujii : Quantum Optical Construction of Generalized Pauli and Walsh-Hadamard 



Matrices in Three Level Systems, quant-ph/0309132, 

[3] K. Fujii, K. Higashida, R. Kato and Y. Wada : Cavity QED and Quantum Com- 
putation in the Weak Coupling Regime, J. Opt. B: Quantum and Semiclass. Opt, 
6(2004) 502, |quant-ph/0407014| 



[4] K. Fujii, K. Higashida, R. Kato and Y. Wada : Cavity QED and Quantum Compu- 
tation in the Weak Coupling Regime II : Complete Construction of the Controlled- 
Controlled NOT Gate, to appear in the book "Trends in Quantum Computing 
Research", 2006, Nova Science Publishers, Inc (USA), |quant-ph/0501046 



[5] K. Fujii, K. Higashida, R. Kato and Y. Wada : N Level System with RWA and 
Analytical Solutions Revisited, quant-ph/0307066 



[6] K. Fujii and H. Oike : How to Calculate the Exponential of Matrices, 
|quant-ph/0406115| 

[7] K. Fujii : A Modern Introduction to Cardano and Ferrari Formulas in the Algebraic 
Equations, quant-ph/0311102 



10 



